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Abstract: Adding matter of mass m, in the fundamental representation of SU(N), to 
N = 4 supersymmetric Yang-Mills theory, we study "generalized quarkonium" containing a 
(s)quark, an anti(s)quark, and J massless (or very light) adjoint particles. At large 't Hooft 
coupling A > 1, the states of spin < 1 are surprisingly light (Kruczenski et al., hep- 
th/0304032) and small (hep-th/0312071) with a J-independent size of order yX/m. This 
"trapping" of adjoint matter in a region small compared with its Compton wavelength and 
compared to any confinement scale in the theory is an unfamiliar phenomenon, as it does 
not occur at small A. We explore adjoint trapping further by considering the limit of large J. 
In particular, for J S> y/\ 3> 1, we expect the trapping phenomenon to become unstable. 
Using Wilson loop methods, we show that a sharp transition, in which the generalized 
quarkonium states become unbound (for massless adjoints) occurs at J ~ 0.22\/X. If the 
adjoint scalars of N = 4 are massive and the theory is confining (as, for instance, in N = 1* 
theories) then the transition becomes a cross-over, across which the size of the states changes 
rapidly from ~ \f\~/m to something of order the confinement scale ~ A -1 . 
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1. Introduction 

In order to bring AdS/CFT techniques to bear even remotely upon QCD, the original 
proposal jl] must be supplemented by somehow introducing quarks in the fundamental 
representation. This was first done || in the system of D3 branes at an orientifold 7- 
plane, which requires the presence of four D7-branes; the field theory is N = 2 Sp(N) 
with four hypermultiplets in the fundamental representation and one in the antisymmetric 
representation. However, the technical difficulties, both in the field theory and in the 
supergravity, of this system, and prejudices about what to use it for, obstructed progress 
for some time. In || Karch and Katz cut through these barriers by pointing out that one 
could in the AdS context simply add a finite number Nf of M = 2 hypermultiplets to SU (N) 
J\f = 4 Yang-Mills, since as N —* oo the positive beta function which results is negligible. 
The details of the gauge theory (and the notation we will use to describe it) are presented 
in Sec. 2 of Ref. [Q]; we will not repeat them here. On the AdS side, this corresponds to 
adding a finite number of D7-branes into the AdS^ x S 5 geometry, and observing that the 
backreaction of the 7-branes on the geometry and the dilaton is a subleading effect. 

This simple observation then led to a number of interesting developments. First, Karch, 
Katz and Weiner || showed that the system exhibits what one might call 'Gribov confine- 
ment' or simply 'strong-field confinement', which involves confinement of heavy quarks 
without the presence of flux tubes. In M = 4 Yang-Mills plus massive matter in the fun- 
damental representation, the absence of flux tubes is evident since the theory is conformal 
in the infrared. However, as shown in ||, if a heavy (s)quark and anti(s)quark of mass 
M are pair-created and gradually separated from one another, and if the theory contains 
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(s)quarks of mass m < M, then under some circumstances a pair of the lighter (s)quarks 
will be created, confining the heavy quarks; and this happens even though the theory is 
infrared-conformal and does not confine electric flux. 

Next, within the same theories, the dynamics of quarkonium states (and generalized 
quarkonium containing additional adjoint matter) were studied in Q. In particular, bound 
states containing a squark Q of mass m, an antisquark Q of mass m, and some number J 
of massless fields <1> from the N = 4 sector were studied, along with their superpartners. 
The spectrum of states for various spins S was studied, with exact results obtained for 
S < 1, which are states described purely by the eight-dimensional gauge theory living on 
the worldvolume of the D7-brane embedded in AdS§ x S 5 . The low-spin states were seen 
to be surprisingly deeply bound, with masses of order m./, = mj\f\ , where A = 4irg s N ^ 1 
is the 't Hooft coupling. 

To explore these deeply-bound states further, we recently @ studied their form factors 
with respect to conserved U(Nf) flavor currents, SO{4) R-symmetry currents, and the 
energy-momentum tensor. To regulate certain computations, we considered adding masses 
m<j> to the $ particles in some cases. Among other observations, we discovered that these 
states all have size of order m^ 1 = ^/Xm" 1 , independent of J and of their radial excitation 
number. Moreover, we found that the states have a size which is not sensitive to the mass 
m<j> of the $ particles (as long as m$ <C mn-) This is very different from the weak coupling 
regime. A QQ state has the physics of the hydrogen atom (and size (Am) -1 .) However, the 
Q&Q state is dynamically more similar to the hydrogen molecule and has a much larger 
size; the length scale of the wave function for Q and Q is a geometric combination of m _1 
and . In particular, at weak coupling, the size of the Q&Q state diverges as m$ goes to 
zero, 1 whereas at large A we found J^j that it remains finite and of order mZ . 

In short, there is a new phenomenon at large A not previously observed in gauge theory, 
in which light particles $ are trapped in a region which is small compared both to their 
Compton wavelength m^ 1 and to the distance scale at which electric flux is confined A -1 . 
This "trapping" effect, which we believe is a new phenomenon and which is related to 
other unfamiliar stringy effects at large A, is what we will seek to explore further in this 
paper. We will argue that when the number of adjoints in the generalized quarkonium state 
becomes parametrically of order the state in the string theory ceases to be a pointlike 
gauge boson on the D7-brane; instead it becomes a long semiclassical string that hangs 
down below the D7 brane. This in turn means that the generalized quarkonium state in 
the field theory is becoming large, and the trapping is becoming ineffective. If the adjoints 
are massless, complete untrapping occurs when the number of adjoints is of order 0.22\/A- 
We will see that this untrapping transition occurs very rapidly as a function of either A or 
J; for J > J* the states are unbound, while for </<«/* the states of low spin are generally 
trapped. If the adjoints are massive and/or the theory is confining, then the untrapping 

1 This can be seen as follows. The Q<&Q state at large N has an attractive Coulomb potential between <!> 
and Q, and also between $ and Q, but no potential between Q and Q (more precisely, there is a repulsive 
but 1/iV-suppressed potential.) A straightforward Born-Oppenheimer calculation shows the Q and Q move 
independently and slowly in a wide and rather flat potential- well induced by the rapidly moving it can 
easily be checked that the overall size of this well grows as me decreases. 
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transition is interrupted when the size of the generalized quarkonium state becomes of order 
the confinement length or the Compton wavelength of the adjoints, whichever is smaller. 2 

Our methods for establishing these claims will be as follows. First, we will argue, both 
on the basis of the spectrum computed in and using a BMN-type argument, that the 
generalized quarkonium states at large J should become unbound at some J = J* where J* 
is parametrically of order The BMN argument confirms, moreover, that these states 
are not metastable for J ~ vA, and that a loss of efficient trapping occurs as J approaches 
vA. At some point this allows the generalized quarkonium states to become much larger 
than their inverse masses. We therefore argue that for J just below J*, the states should 
be relatively large in size and have small binding energy. Such states would resemble a 
hydrogen molecule, with fast light adjoints weakly bound to two heavy (s)quarks. We 
therefore expect the motion of the (s)quarks to be non-relativistic, for J sufficiently close 
to J*, and therefore a Born-Oppenheimer-type calculation is appropriate, in which the fast 
motion of the adjoints is treated first, allowing for the computation of an effective potential 
in which the heavy (s)quarks move, and the slow motion of the heavy (s)quark states in this 
potential is treated second. The computation of the effective potential is simply a Wilson 
loop computation, which we carry out in section 3. This calculation has its own interest and 
we explore it in some detail. As expected, the Wilson loop shows a gradually decreasing 
level of trapping, with trapping entirely lost at J = J* = 0.22\/X The result allows us 
to conclude that the effective potential in the Born-Oppenheimer calculation for the states 
with dynamical (s)quarks has the same transition; above a critical value of J, the effective 
potential for the (s)quarks is simply zero, and no generalized quarkonium bound states can 
form. Just below this value of J, the effective potential for the (s)quarks is Coulombic with 
a small coefficient, so the orbits of the (s)quarks are hydrogenic, with a computable effective 
coupling. This establishes that our picture for the loss of trapping and the unbinding of the 
states is self-consistent. 

2. Trapping Many Adjoints? 

The fact that the trapping effect creates generalized quarkonium states Q& J Q whose size 
does not grow with J raises a question as to what happens when J becomes large compared 
to \/A- Supergravity cannot describe these states; instead they are expected to be better 
described using a string theory obtained in a BMN limit Q|. Most work on BMN limits has 
attempted to describe operators in conformal field theories (or states of field theories on 
spatial three-spheres) but a BMN limit describing massive states of nonconformal confining 
field theories in four-dimensional Minkowski space was obtained in ||. Could a similar 
BMN limit describe the massive S = 0, f states of the form Q& J Q for large J? 

There is a simple reason to suspect the answer is no. According to [0, the ground state 
consisting of one Q, one Q and J massless $ particles has mass of order Jrrih = (J /\f~\)m. 
Clearly when this is greater than 2m, the system should not be bound. However, this 
reasoning needs to be checked, especially as large 't Hooft coupling physics has often held 

2 Our conclusions regarding the states with large numbers of adjoints differ significantly from the prelim- 
inary suggestions made in M. 
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surprises. In particular, the argument could be correct, but still there could be a potential 
barrier which makes these states metastable and forces them to decay via tunneling. To 
explore this latter possibility we have computed the BMN limit corresponding to these 
states. 3 As we will see, there is no metastability. 

Our setup is that iV D3-branes fill the 0123 directions of the ten-dimensional space, 
and are located at the origin of the 456789 coordinates. A D7-brane probe is placed at the 
position x 8 = vrtQCt = L, x 9 = 0, and fills the 01234567 directions. We can write down 
the metric in a form that manifestly shows the embedding of the induced metric on the 
D7-brane 



ds 2 = a' 



p2 + p2± v ^dx»dx» + -^^{dp 2 + P 2 dn 2 + dp 2 ± + p\d x 2 



(2.1) 



v/A ^ ' p 2 + p\ 

where x^, p = 0,1,2,3 are the ordinary field-theory Minkowski space coordinates, and 
the remaining part of the D7-brane world-volume is spanned by X ■ X ■ X } X , which we 
write using rescaled spherical coordinates p 2 = (X)i=4 X T)I a ' 2 an d d£l 2 = dip 2 cos 2 9 + d9 2 + 
sin 2 9d£l\. The space transverse to the D7-brane is represented by rescaled polar coordinates 
p\ = [x\ + Xg)/a' 2 and x- I n these coordinates, p and p± have mass dimension +1, while 
Minkowski space coordinates x^ have mass dimension —1 as usual. 

To find a BMN limit for states of definite mass (static states which are eigenstates of 
the field-theory Hamiltonian d/dt) we should seek to take a Penrose limit with respect to 
a null geodesic at a constant AdS radius, which we should expect to lie near p = 0. This is 
because H hadrons of high charge correspond to modes which are concentrated close to a 
nonzero and small AdS radius. 4 

We therefore seek a null geodesic at a constant point in physical space x 1 ,^ 2 ,^ 3 and 
moving in time and around an equator of the S 3 , e.g. the curve 9 = 0. A particle moving 
on such a curve will have both large energy E and large charge J. The effective Lagrangian 
for a particle describing this kind of motion is 



C = a' 



2 



P 2 + P±12 , VX 2 . 2 V\p - 2 



p 2 + p\ p 2 + pI 

where the • means the derivative with respect to the affine parameter. The null condition 
C = gives 

where E and J are the conserved energy and angular momentum associated with the Killing 
vectors d/dt and d/dtp respectively. This is the dynamics of a particle moving in a potential 
V(PiP±) = P + p\l P- Note that the D7-brane is located at p± = uiq. If the string itself 
were fixed to be at p± = uiq, then the minimum of the potential would be at p = p± = vtlq. 
However, the string (except for its ends) is free to move to any value of p±, and with no 
constraint on p± the minimum of the potential lies at p± — > 0. This shows that there is 
no stable geodesic with the properties that we are seeking. In particular, a particle, or an 
unconstrained piece of string, on the above trajectory wants to fall off the D7-brane toward 
the horizon of the AdS space. 

3 This computation was also discussed briefly in 

4 A hadron's wavefunction falls off as r~ A , where r is the AdS radius and A is the dimension of the 

lowest-dimension operator which can create the hadron. A hadron of large charge J can be created only by 

an operator of large charge which, since it contains of order J fields, has A ~ J. Therefore, the wavefunction 

for a hadron of large charge has a narrow peak at a small AdS radius. 
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Figure 1: The proposed transition 
as the R-charge J (the angular mo- 
mentum in the horizontal plane) is 
taken large. The D7 brane, shown 
schematically, is warped by the ge- 
ometry in the radial direction; grav- 
ity pulls down in the figure, along the 
coordinate U = a//? 2 + p\ , and the 
AdS horizon is shown as a horizon- 
tal plane, (a) A small string with 
moderate J orbits around the end of 
the D7 brane. (b) As J — > J* the 
string distends and moves off the or- 
bit shown in (a), (c) The ends of 
the string descend to the origin of 
the D7 brane, while its middle ex- 
tends below the D7 brane and carries 
the angular momentum J. (d) At 
J = J* the string reaches the hori- 
zon, which absorbs the angular mo- 
mentum, leaving two long stationary 
strings — unbound (s)quarks. 



This result shows the Q& J Q states are fully unsta- 
ble, not metastable, for J 3> vA. In particular, we are 
led to a speculation concerning the nature of that in- 
stability. These states correspond to open strings with 
both ends on the D7-brane; while the ends cannot move 
off the D7-brane, the bulk of the string worldsheet may 
be expected, according to the above computation, to 
fall off the D7-brane toward the horizon. This is illus- 
trated in Fig. |j. Once the string becomes an extended 
object, and no longer acts as a massless point particle 
as viewed in ten dimensions, it will no longer follow the 
orbit shown in Fig. ||a. (This represents the onset of 
the loss of efficient trapping in the gauge theory.) In- 
stead, the string will move to a lower orbit, with more 
of the angular momentum being carried by the lower 
regions of the string where the energy cost is less due 
to the redshift of the metric. For large enough J, the 
ends of the string will move almost to the origin of the 
D7 brane, with the angular momentum mainly carried 
by the center of the string, which continues to orbit 
around the origin. As J increases further, the string 
will become longer still, and eventually the string will 
hit the horizon, at which point the angular momentum 
will dissolve into the horizon, as shown in Fig. |l]d. In 
this case the two (s)quarks, represented by two indepen- 
dent straight strings, are unbound, and the J adjoints 
are free to move away from the (s)quarks and from each 
other. 5 

For just slightly smaller J, corresponding to the 
limit of Fig. 0c just before it transitions to Fig. |l]d, 
the binding of the generalized quarkonium state is very 
weak. In this regime a Born-Oppenheimer calculation 
will be valid, in which the ends of the string are held 



fixed, a Wilson-loop computation is performed first, and 
the motion of the ends of the string is quantized second. This then motivates us to carry 
out the Wilson-loop calculation, which should show the unbinding transition at J = J*. 

3. A Wilson Loop and the Unbinding of Adjoints 

The Wilson-loop computation in field theory corresponds, as is well known, to the problem 
of computing the energy of an appropriate semiclassical string sitting within a supergravity 
background. We proceed in perfect analogy with the hydrogen molecule calculation of the 
effective potential between the two p rotons induced by the electrons. Instead of studying 
5 To be precise, the colors of the (s)quark and anti(s)quark are uncorrelated, because they interact at the 
planar level only through the adjoints; as the adjoints move away, the residual force between the (s)quarks 
is only of order 1/iV 2 , and is repulsive. Note also that the Minkowski spatial dimensions are of necessity 
suppressed in Fig. the (s)quarks are not at the same point in space, and the force between them is finite. 
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(s) quarks of finite mass m at a distance of order 1/m, we will take m to infinity while 
fixing the distance L between the heavy squark and antisquark. Meanwhile we keep m$ 
fixed. If we are in the regime where the $ particles are trapped, and L < mj , then 
the $ particles should spread out over a region whose size is of order L. If they are not 
trapped, they should spread out over a region whose size is of order m^ 1 . We will see a 
transition between these two behaviors below. The results of this calculation will give us 
the effective potential. This potential will remain valid for (s)quarks of finite mass m if 
the binding energy Eq of the ground state for finite-mass quarks in this potential is small 
compared to m. For m<j> sufficiently small, the potential must be of Coulomb type (by 
conformal invariance). A necessary condition for small binding energy is that the coefficient 
of the Coulomb term — the effective a — be small compared to 1. This is always true 
sufficiently near the transition point J = J*, and this will allow us to compute J*. Note 
however that for any fixed A there need not exist any bound states well-described by the 
Born-Oppenheimer method, and indeed we will see that such states are not generic. 

The easiest way to study this problem initially is to take m$ to zero. In this case we 
need only make appropriate modifications of the original computation of the Wilson loop 
by Maldacena |l0| and by Rey and Yee [11| to determine the coefficient of the 1/L potential 
between infinitely heavy sources in the fundamental representation of SU(N) in J\f = 4 
supersymmetric Yang-Mills. It was found that the coefficient was -4vr 2 VA/[r(l/4)] 4 when 
the sources are antiparticles of one another. Maldacena also considered the case when the 
sources couple differently to the scalar fields of the N = 4 theory, a possibility we will return 
to in just a moment. 

To understand the calculation we wish to perform, we need to examine the theory 
carefully. The N = 4 theory has an SO (6) R-symmetry, and a sextuplet of scalar fields 
which it is convenient to organize into three complex scalars <3?j. To preserve the maximal 
supersymmetry, a BPS-saturated source in the fundamental representation of SU(N) must 
be a massive vector multiplet, which preserves an SO (5) subgroup of this symmetry. Al- 
ternatively, we can preserve half the supersymmetry if we add a source which is a massive 
J\f = 2 hypermultiplet; this source preserves an SO (4) subgroup of the R-symmetry. In 
both cases, the source-antisource state appropriate to a Wilson loop computation breaks 
supersymmetry (in general), but preserves the same SO (5) [or SO (4)] subgroup as the iso- 
lated BPS source. In the computation of [10, If |, this source-antisource state is described as 
a string whose ends lie at the boundary of AdS a distance L apart, and which are oriented 
on the 5 5 so as to preserve the appropriate SO (5) symmetry [or SO (4).] 

We may now ask that the bulk of the string worldsheet be allowed to rotate around 
an equator of the 5 5 , so that it picks up a charge J with respect to an 50(2) « U(l) 
subgroup of the 50(5) [or 50(4).] This corresponds to the source-antisource state binding 
to J complex scalars <3? which are charged with respect to this U(l). (Calculations of this 
sort of state have also been undertaken in various papers, especially in |I3|, |14|, where a 
very similar method was needed.) For the case relevant to the D3-D7 system we have been 
considering in this paper, we imagine we introduce an infinitely massive hypermultiplet 
Q, Q which couples to the scalar field $3 (using a D7-brane at x 8 = m — > 00, x 9 = 0.) We 
may then add J scalars <3?i to the state built from a (s)quark and anti(s)quark by allowing 
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the string to rotate in the x 4 -x 5 plane. 

Since the theory with m$ = and m = oo is conformal N = 4 Yang-Mills, we know the 
potential will be of the form V(L) = /(A, J)/L. Our simplest goal is to determine /(A, J). 
More generally we wish to determine the shape of the string corresponding to a given A and 
J. We will now present this calculation, which exhibits a trapping transition. 

The computation is slightly more complicated than the prototype [[U], O, since a third 
coordinate comes into play. However, the mathematics largely reduces to one of Maldacena's 
other computations, as we will see. A even more similar computation was performed by 



Tseytlin and Zarembo in [14|, and our techniques follow theirs. On the surface it would 
appear that we are about to repeat their computation, but the details of their solutions are 
crucially different. In particular, we choose different boundary conditions from [|l4"|| . Our 
boundary condition corresponds (as we will see) to a distribution of the global U(l) charge 
which is regular near the source and antisource; that of |l4j is singular (though integrable.) 
The instability found at weak-coupling in [14| does not apply to our computation. 6 

The equations for the shape of the string and its energy are mathematically equivalent 



to the computation of [10] in which the source and antisource are not each other's antiparti- 
cles. Maldacena introduced two massive vector multiplets which preserved different SO (5) 
symmetries, leaving an SO(4). An angle Ax enters this computation, describing how badly 
aligned are the SO(5)s preserved by the source and antisource; more intuitively, it specifies 
the angular separation on the 5 5 of the two ends of the string as they approach the bound- 
ary. For Ax = the source and antisource are antiparticles, while for Ax = vr the state 
of the source and antisource is BPS-saturated and the binding energy is zero. In this case 
V(L) = /(A, Ax)/L; Maldacena found an implicit form for the function /. We will see that 
this function reappears in the calculation below, but with a very different interpretation. 

3.1 The Calculation 

We will consider a rectangular Wilson loop of length L and duration T, and we consider the 
limit T — > oo in order to extract the potential energy between sources at a distance L from 
one another. This calculation is dual to a computation of the energy of a string whose ends 
lie on the boundary of AdS and are separated by a distance L in the spatial coordinates of 
the gauge theory. 

We can use the Polyakov action to describe the semiclassical string worldsheet 
S = ~4^b / drdai-^/^GMNd^d^ 

with appropriate boundary conditions at the end points, which lie at p± = oo. 

We are looking for a stationary configuration of the string such that x = 0- The string 
should be rotating around an equator of S* 3 , parameterized by the angle ipj we set ip = ut 

6 The Euclidean-space calculation of solution (B) in Sec. 3.2.2 of jl4| has mathematical similarities to our 
solutions, but arises from a different boundary condition (again containing a singular charge distribution) 
and has a correspondingly different interpretation. To be more precise, the computations in G^] are done 
using a boundary condition that a coordinate called "t/>" goes to at the boundary. Our coordinate 8 is 
■| minus this u ip" , and we instead choose the condition 9 — > at the boundary; this minimizes the global 
charge density near the source and antisource and reduces the energy. 
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where uj is a constant. We put the (s)quark (one end of the string) at X(p± = oo) = L/2 and 
the anti(s)quark (the other end) at X(p± = oo) = —L/2, in order that string configuration 
be symmetric about X = 0. We will find it useful to employ the coordinates U and 9, 

U 2 = p 2 + p 2 < , 9 = arctan — . 

P± 

U is the usual AdS§ radial coordinate, and 9 is a polar angle on the five-sphere. We are 
interested in solutions in which the source itself is not associated with any global charge, 8 
and so we expect 9 that 9 — ► and X(U) goes to the solution of [10, [Tl]] as U — > oo. 



Using the conformal gauge 7 Q/ g = r] a a , the required portion of the string action (setting 
the string tension equal to 1) becomes 

T f 

S = — — / daC = TSq 
4vr J 

c = p -^{? + x' 2 ) + -^(p 12 + p'I - pH 2 ) , (3.1) 

VA P 2 + Pi 

where ' = d T and ' = d a . Using the coordinates U and 9, we can separate C into two parts 
which depend only on the AdS§ coordinates and the S 5 coordinates, respectively. Fixing 

T = t, 

C = ^(1 + X' 2 ) + ^U' 2 + V\(9' 2 - lo 2 sin 2 9) . (3.2) 

The equations of motion for the AdS§ coordinates U, X decouple from that of the S 5 
coordinate 9. 

X' = gU 2 /U 2 (3.3) 

rr/2 2 tjA i 2 _ i 

U | 9 U Q 1 rr2 _ 9 l _rj2 ( oa\ 

9' 2 + w 2 sin 2 9 = uj 2 sin 2 O • (3.5) 

Here g is a (dimensionless) constant of integration. At the midpoint of the string X = 0, 
where 9' = and {/' = 0, 9 reaches its maximum value 9q and U reaches its minimum value 
Uq. The constraint equation resulting from the conformal freedom of the worldsheet metric 
in the Polyakov action, 

^(-1 + x' 2 ) + ^U' 2 + V\(9' 2 + u 2 sin 2 9) = , 



7 With the coordinate change, the metric (2.1) can be written as 

J 2 / 

as — a 



rj2 rr 

^i-dt 2 + dx 2 ) + i-rdU 2 + V\(de 2 + sin 2 6 d<b 2 + cos 2 9 dy 2 ) 
s/\ U 2 



which shows that 9 is a polar angle on the full S 5 . 

8 At most, the global charge of the source should be of order 1, not of order J, in order to correspond to 
an infinitely massive hypermultiplet. 

9 Note the subtleties addressed in 11 211 do not arise here. 
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gives the relation 

2 



(f - 1 - ft 2 sin 2 9 , n = . (3.6) 



The state we are studying has angular momentum J and energy E, 

J=^°=^ [da sin 2 9 (3.7) 

E = u^-S = [ daU 2 . (3.8) 



5uj 2tt y/\ , 

Note that the energy is divergent and must be regulated; we want the negative potential 
energy between the source and antisource, so we must carefully subtract the divergent 
masses of the sources, as in 11 1. 

Remarkably, part of the solution to these equations is of the same form as one of 
Maldacena's computations. In particular, the configuration X = X(U) is the same as 
that of a string stretched between two D3-branes located at the boundary with an angular 
separation on the S 5 if we identify g 2 in our equations with 1 — I 2 in equations (4.10)-(4.12) 
of [ 10 1 . The string configuration 10 X = X(U) is obtained from Eq. ( |3.4[) 

X = ^-I(g,U/U ) , (3.9) 

where 

rU/Vo i 

I(9,U/U )= dy (3.10) 

Ji rv(y - i)(y +g ) 

= 1 [E(wr\U /U),ig)-F(Bm-\U /U),ig)] . 

(Here E, F and later K are elliptic integrals, with conventions defined in Appendix C.) The 
string end point is located at (X, U) = (L/2, oo), which determines Uo for a given g: 

where 

1(g) = I(g, U/U - oo) = \ \E(ig) - K(ig)\ . (3.12) 

T 

The relation between 9(U) in our computation and the mathematically-related angular- 
separation variable in |l(| is not so direct, however. The string configuration 6 = 0(U) is 
obtained from Eqs. (|3.5|) and (|3.9|): 



./l _ Ju/u \/{y 2 - l){y 2 + 9 



3 Since the string configuration is symmetric about X = 0, we describe only the X > half of the string. 
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Here we used 9{U — > oo) — ► 0, consistent with the boundary condition. 11 While the right- 
hand side of this equation appears verbatim in equation (4.10) of [pi], the left-hand side is 
significantly different in all respects. 

Equation ( |3.13j ) may also be written as 



sin0 o F^in-^sinfl/smflo^smfloJ = V 1 ~ 9 2 F (sin- 1 (U /U),ig) . 

The condition that 9 = 9q at U = Uq gives the relation which determines g for a given 9q: 

sin 6»o K(sm9 ) = y / l-g 2 K(ig) . (3.14) 

Using the equations of motion, the angular momentum J can be written as 

,2 a 



k Jo 



sin 



TT 



sin f — sm 
K(sm9 ) -E(sm9 



(3.15) 



The energy in Eq. ( |3.8| ) should be regularized by subtracting the masses of Q and Q because 
it includes infinite Q and Q masses. 



EL 



TT 



dy 



y 



V(y 2 -i)(y 2 +5 2 ) 



1 



2VXg 3 



TT 



(3.16) 



which again matches equation (4.13) of [10|, but with quite a different interpretation, as we 
will now see. 

Our solution is now complete: we can make one-to-one correspondence between the 
parameters in our equations with the physical quantities, J and E. For a given J, there 
is a corresponding 9$ from the relation Eq. ( |3.15| ) and g is determined from Eq. ( p. 14 ). 
Accordingly, E is obtained from Eq. ( 3.16[ ). Finally, Uq is determined from the relation 
Eq. (3.11), and U(X) can be implicitly found from Eq. (|3.S|), 

However, it can be seen that the range of J is restricted, and cannot be arbitrarily 
large. This is the sign of the instability which we have been seeking. The left side of 



Eq. ( 3.14 ) is an increasing function of 9q, while the right side is a decreasing function of g. 
The right-hand side of Eq. ( |3.14 ) has an upper-bound of vr/2, which occurs when g — > 0, 
while the left-hand side diverges logarithmically as #0 ~~ > vr/2. Therefore, there exists an 
upper bound for 9r>, which we will call 0*, with 



sin#* K (sin 9*) 



TT 



(3.17) 



11 Curiously, the equations also allow for a solution with 9 — ► 60 as U — > 00, one which has the same 
relation between g and Oo and which has the same J and U(X). However, the energy of such a state is 
infinitely larger than the one of interest to us, reflecting the fact that a source with fixed nonzero 9 is a much 
longer string than one with 9 — fixed. 
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Figure 2: String configurations for g — l(dotted), 0.1(thin solid), and 0(thick solid) (equivalently, 
J I J* = 0, 0.9867, and 1, respectively) in the X-U plane. Note the minimum value of U is a 
decreasing function of J. 



and correspondingly one for J, which we will call J* 



7T L 



K(sin0*) -E(sin6* 



Numerically, sin6»* ps 0.793 and J* ps 0.22\/X . 

To know what does happen when J — ► J*, we need to analyze the behavior of the 
string 12 as g — > 0. Since /(g) — ► 7r/4 as g — > (J — ► J*), 



-5 



. 



(3.18) 



This shows that Uq touches the horizon and the interaction energy vanishes as J —* J*. 
That is, for fixed L, i/ie massless $ particles completely unbind from the infinitely heavy 
(s)quark and anti(s)quark for J > J*. As we will discuss further below, this implies that 
generalized quarkonium states also become unbound (for massless adjoints) for J greater 
than J*. 

The relation between E and J (for fixed L) is plotted numerically in Figure || and is 
obtained in two limits ( J — > and J — » J*) in Appendix A; the result is 



EL 



EL 



4tt 2 V\ 

[r(i)] 4 



-n I - 



3/2 



J 

2tt^ + 



Jin J 

y/X sin 9* 



J^O 



3/2 



+ 



(3.19) 
(3.20) 



Note \E\ decreases as J increases, as expected; the presence of the scalars in the state 
makes the binding energy smaller in magnitude. Interestingly, the relation between E and 
J is nonanalytic as J —* J*, indicating a sharp transition. 



12 In this limit, uj — > (7o/(v / Asin6'*) ~ g — » 0. Note w is not a physical quantity; it is J that is physical, 
and larger J does not necessarily means larger w. 
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EL/ VI 




It is interesting to ask where the scalar fields actually lie in the space between the source 
and antisource. We do not know the precise answer to this question, but some significant 
insight is obtained by examining the density of J as a function of x = X/[L/2], or as a 
function of y = U/Uq. For a given g, the density of J is given by 



p{y) « 



sin 2 9{y) 



V(y 2 -i)(y 2 + 9 2 ) 



p(x) oc y 2 sin 2 %) , 



(3.21) 



where 6{y) is given in Eq. ( |3.13 ). This is illustrated in Figure |j. As can be seen, as J — > J* 
the scalars cluster at U ~ Uq — > 0; using the usual inverse relationship between the U 
coordinate and the size of distributions in the field theory [jl], ft5j, this strongly suggests 
that the scalars spread out, and thus the size of the state is diverging, as J — > J*. 13 



p[y]/J 



pM/J 

0.7~ 



(a) 



(b) 



Figure 4: (a) Density of J versus y = U/Uq and (b) Density of J versus x — X/(L/2) inEq. ( ^2l| ) 
for g = 0.1 (J/J* w 0.9867.) 



In summary, we have found that in the case of infinitely massive (s) quarks and massless 
$ adjoints, trapping occurs for J < J* 0.22\/A, decreasing in its effects as J — > J*. For 
J > J*, the adjoints do not bind to the sources, and the effective potential vanishes. 



13 In the Minkowski-space solution of [E3, the factors of sin 2 9{y) are replaced with 1. This leads to more 
singular distributions of charge, since sin (9 ~ 1/y at large y. 
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We have also considered how this calculation is modified if, as in [10|, the string's 
endpoints are at two different positions on the five-sphere. Without loss of generality we 
may take these positions to lie on the same great circle, one of them at +Ax/2 and one at 
—Ax/2. This computation, presented in Appendix B, shows qualitatively similar behavior; 
as A% increases, the value of J at which the binding energy drops to zero decreases, until 
at Ax = 7r, where the state is BPS saturated, the binding energy is zero even for J = 0. 
The curve in the (J, Ax) plane where the binding energy vanishes is shown in Figure ^. 
Interestingly, the U (X) curve and the binding energy are functions of only one combination 
of J and Ax, given by the parameter g; in particular, for every J < J* with Ax = 0, there 
exists a value of Ax with J = that leads to the same curve U (X) and the same binding 
energy. We do not know of a deep reason for this. 



3.2 Discussion 

Let us now consider what the Wilson loop computation implies for the generalized quarko- 
nium states, for which m$ is not necessarily zero and m = vtiq is finite but much larger 
than m$. It is useful to break the discussion up using a couple of intermediate steps. 

For infinite m but nonzero m$, we expect the new mass-scale to introduce some new 
physics at a radial position U oc m<j> in the AdS space. Since the theory is no longer 
conformal, the potential will now be V{L) = /(A, J; m$L)/L. For instance, if all six scalars 
receive positive masses, as in the Af = 1* theory fl~6|| , then we expect the AdS space is 
effectively cut off at U m in ~ m®- This means that our solutions will be valid only until 
Uo reaches m<j>. This in turn means that the shape of the string, and the corresponding 
effective potential, will change once L > mj . Indeed, the calculation should eventually 
match on to that of ||. For large L, the string will lie at U = U m i n for most of the region 
—L/2 < X < L/2, and we expect the ^-coordinate of the string to relax to 7r/2, as this 
should minimize the energy of the configuration. This gives || 

E « ^T 2 L 2 + m|j 2 (U min ~ m*) . (3.22) 

Here T ~ 1S the tension of a standard flux tube in the gauge theory. Thus, the 

potential becomes linear in L for sufficiently large L and fixed J, and becomes linear in J 
for fixed L and sufficiently large J S> VX The regions with different behaviors are shown 
in Figure [|. Note all the transitions are cross-overs; there are no phase transitions once m<i> 
is nonzero. 

For finite niQ and m$ = 0, on the other hand, the static Wilson loop we have been 
studying is no longer physical; instead the (s)quark and anti(s)quark must be put in orbit 
around each other. We have not attempted to examine this more complex dynamical 
problem carefully for general J, but in the regime J ~ J* the situation is under control. As 
we have discussed earlier, in this regime the force between the (s)quark and anti(s)quark is 
small and their motion is slow and nonrelativistic, justifying a Born-Oppenheimer approach. 
The effective potential for the slow (s)quark and anti(s)quark is given by the Wilson loop 
calculation we just performed. As long as a e ff = EL, as given in Eq. (EOSD , is small 
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j/Va 



0.6 



0.4 



0.2 



rruL 



Figure 5: The various regimes for nonzero m$. The diagonal dashed line is J = ^m<$,L. The 
potential goes as 1/L in the lower-left region, as Jm$ in the upper left region, and as L times 
the flux-tube tension (of order ^/\m%) in the lower-right region. All borders between regions are 
cross-overs; there are no sharp transitions. 



compared to 1 (as becomes true in the regime J « J*), then we would expect that the 
heavy particles will move with a velocity v ~ a e ff <C 1, as required. 14 

However, to see whether and when this computation is valid, we must self-consistently 
solve for L, given a fixed quark mass mg. The scaling of the various quantities as g — > 
is now different from the scaling relation ( |3.18D , which was obtained holding L fixed. Since 
L = (a e fjmQ) _1 in a Coulomb system and Uq ~ VXg/L, the scaling of the various quantities 
with g ~ (J* - J) 1 /^- 1 / 4 as g -»■ is 15 

a< s = EL~\ 1 /*tf i , L-A-V, E~\g\ U ~ Xg 4 . 

If we only required a c s <C 1, we would allow J* — J <C A 1 / 4 . But this is not sufficient. We 
must also require Uq <C wq, so that our Wilson line computation with strings that extend 
to infinity bears some relation to the calculation with strings that extend only to the D7 
brane at U = vtlq. Furthermore, we must also require L S> yX/mQ, the trapping scale, 
which we showed in [|j also sets the scale of the (s)quark and anti(s)quark wave functions 
for the trapped states. The last condition requires a e g <C l/\/X -C 1. This in turn requires 
J* - J ~ A" 1 / 6 < 1! This means that for most values of A there is no integer J for which 
there is a state described by a Born-Oppenheimer calculation. 

Still, this scaling is still sufficient for our purposes: for fixed J, there is always some 
range of A for which the calculation is valid. If we follow any particular state of definite 
J as A is decreased adiabatically, the state becomes unbound at J = J* by passing, rather 
swiftly, through a regime in which it is described by our Born-Oppenheimer calculation, 



14 There will be a full Coulomb spectrum of states in this regime. Of these, only a handful are adiabatically 
related to the states with adjoint trapping, in particular only those states with spin < 1 and radial excitation 
n with J + n < A 1 ^ 2 jjj. The rest of the states are never described as supergravity states for any A, requiring 
instead the full string theory; correspondingly they are always much larger in size than the trapped states. 



See Eq. (A..4) in the Appendix A. 
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grows much larger than its typical size y/X/vriQ, and unbinds when J* ~ y/\ reaches J. 
Consequently, for any fixed A, there are no bound states with J > J*. All states with low 
spin, no radial excitation, and J* — J > A -1 / 6 have adjoints trapped in a small region. This 
is consistent with our earlier conjectures; what we have learned with these calculations is 
how rapidly the transition occurs. 

Finally, let us consider what will happen if both m$ and mg are nonzero, with m$ <C 
m,Q. In this case we expect a cross-over, as in the Wilson- loop computation, to occur in 
the regime J ~ J*. As before, the Born-Oppenheimer regime exists at J* — J ~ A -1 / 6 . 
The calculation of the Coulombic potential and its coefficient need only be modified by 
confinement effects at an even smaller value of J* — J, at the point where the untrapping 
of the adjoints allows them to reach the confinement scale A ~ l/m$. 16 Again there is a 
cross-over, a rather sharp one, in which, for fixed J and decreasing A, each bound state 
goes from a object of size ~ 1/rriQ with trapped adjoints to a larger confined state of a size 
~ l/m$. Equivalently, at any fixed A, there is a sharp transition in the spectrum where 
low-spin zero-radial-excitation states with J < J* have the adjoints trapped in a region of 
size ~ l/ m Q> while states with J > J* spread their adjoints over a region of size of order 
l/m$, with the (s)quark wave functions somewhat more compact. 

Of course, all of these results are modified at finite N. We expect the potential energy 
does not change dramatically. A more important effect is that the strings of moderate J 
(1 < J< \/A) and L <C l/m$ can decay, by emission of closed strings carrying nonzero 
charge J', to open strings with charge J— J'. In other words, generalized quarkonium states 
can decay, via emission of states of the N = 4 or N = 1* theory carrying the global U(l) 
charge. The typical closed string carrying charge J' will correspond (in a confining theory — 
recall that the confinement scale A ~ m$ for large A) to a hadron with mass of order J'm$. 
Meanwhile, for a string with endpoints separated by a length L, dE/dJ ~ 1/L > m$, 
implying that these decays are kinematically allowed. The same is true for dynamical 
generalized quarkonium, for which L ~ \f\m^ when J is not near J*. The phase space for 
these decays is substantial, but by dialing the coupling 1/N the widths of the generalized 
quarkonium states can be made arbitrarily small compared to their masses. We would 
therefore expect them to remain as sharp resonances for large but finite N. 

Finally, let us address the issue of how the physics of large A matches on to that at 
small A. Since trapping occurs for J < J* ~ \/A, it simply need not occur for theories 
with A < 1. This agrees with our understanding of these generalized quarkonium states 
in perturbation theory, which suggests that their size should be of order mj 1 and that 
they should become unbound as m$ — > 0. However, it is worth considering the possibility 
that the J = 1 state in a QCD-like theory with additional adjoint matter might, in very 
favorable circumstances, exhibit adjoint trapping. The quark masses rriQ would need to 
lie not far above A (so that A(mg) ~ 1) with the adjoint masses lighter than A, and 

16 More precisely, for the above statements to be true, we require m$ <JC mA~ 5 / 6 . Moreover when this 
condition is satisfied, the calculation's validity extends somewhat further, since the (s)quark wave functions 
are smaller at this transition point, their size being set by a geometric combination of 1/m* and I/itiq. 
However the physics of the states beyond the transition point as we take J > J, is more complicated, and 
we will not discuss them further here. 
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perhaps additional binding interactions (such as a Yukawa interaction between the quarks) 
might also be needed. This possibility could be explored numerically, though rather large 
N might be required in order to stabilize the state against decay. Although a long-shot, 
the observation of adjoint trapping in lattice gauge theory simulations would certainly be 
remarkable. 
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A. The behavior near J = J* 

To get the behavior of J near J* as g — ► 0, we expand each side of Eq. ( 3.14| ) with respect 
to e = sin 9* — sin 6q and g 2 respectively 

(sin^-6) J Sr(sin^-e) = v / r^/ = = = = (A.l) 

Ji V(y w +9) 

where sin 9* is determined by Eq. ( |3.17| ) . After expanding, 

3vr/8 



If (sin 0*) + sin 0* K' (sin 0*) 
Expansion of J near J* gives 



9 z + 0{g*). (A.2) 



J,-J = e — [if' (sin 0*) - £'(sm0*)] + 0(e 2 ) . (A.3) 



7T 



Therefore, as J — > J* (g — > 0), 



37T/8 j 

-©(sin 0*)/ cos 2 0* 



sine/* — sin #o - — _ 2 - g 



sin0*£'(sin0*) , 
J* — J ~ (sin0* — sin 0o J 

7T COS z #* 

-U) ^Ufc^J ■ (A - 4) 

It is straightforward to get the behavior as J — > (5 — > 1), 

si ^o = I ^§ ] ^(l- 5 ) 1/2 + 0[(l- 5 ) 3/2 ] 

J = ^sin 2 o + O(sin 4 o ) 
£L « + C«/ • (A.5) 
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where E L = -47r 2 VA/[r(l/4)] 4 and ( = 8vr 3 /[r(l/4)] 4 
From the expansion of I(g, y) for small g, 



i(g,y) 



V / 2(^T) + Ob 4 , (y- 1) 3/2 ] , 

(A.6) 

lf-3^-¥ + o(5 4 ,i/y 5 ) , y-oo 

where y = U/Uo, we can get the shape of the string near U = Uq or near the boundary 
(U /Uq — > oo) as J — > J*, 

x 2 = f (l - ^) (y - 1) + Ob 4 , (y - I) 3 / 2 ] , y-1 (A.7) 

l-x = ^fl + ^l + 0(/,l/y 5 ) , y^oo (A.8) 



where x = X/(L/2). 



B. String Stretched Between Two D7-Branes at Different Angles 

We consider a stationary string configuration rotating on S 3 and with its ends on two 
parallel D7-branes, located at U = 00 with angular separation Ax (here X is the polar 
coordinate in the x 8 — x 9 plane, as in the metric (|2.1]).) 
After using conformal gauge and fixing r = t 

C = ^=(1 + x' 2 ) + ^U' 2 + V\(9' 2 + cos 2 9 X ' 2 - co 2 sin 2 9) , (B.l) 
v A U 

where ' = d a . Again the equations of motion for U,X, which are the same as for Ax = 0, 
decouple from the equations of motion for 9, x- 

X' = gU 2 /U 2 (B.2) 
U' 2 ,9 2 U* I 2 _g 2 -l 2 

W + JTp-\ u - — u ° (R3) 

, M cos9 1 

x = ^ (R4) 

9' 2 + lo 2 sin 2 9 + I 2 U ° ™f 9 ° \ = u 2 sin 2 9 + l 2 ^ , (B.5) 
A cos^ 9 A 

where 5 and I are (dimensionless) constants of integration. #0 and Uq are the values of 9 
and U at the midpoint of the string, X = 0, where we take x = 0. At the endpoints of the 
string X -» ±L/2, -» 0, -» 00, and x ±Ax/2. 
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The constraint equation resulting from the conformal freedom of the worldsheet metric 
in the Polyakov action is 



El 



(_1 + a/2) + j^U' 2 + V\(9' 2 + cos 2 6 X ' 2 + ^ sin 2 , 



This gives the relation 



g z + l l = 1 - n 2 sin 2 







U 



As before, the state has angular momentum J and (divergent) energy E, 

5Sq uj^/X 



J 



2ir 



da sin 2 



5S 1 1 f 2 

E = UJ lJ- So = ^7 = xJ daU 

The string configuration is described by the equations 17 : 
cos 9 



d9- 

" • /(I 



-)(eos 2 + 



y/l-l 2 -9 2 



sin 2 O 7 vv ~° " ' ^ 2 



tVt/o \/(2/ 2 - !)(?/ 2 +5 2 ) 



which can be written as 



v^ 2 + Z 2 



sin 0o F sin 



-i 



8111 ( 



sin uq 



O sin 0o 

' vn 2 + z 2 



(B.6) 



(B.7) 
(B.8) 

(B.9) 



^/l-g2-p F ( sin -l(u /U),ig) , 



Ax 



X 



Z cos 0o 

nsin9 J cos 0/ (1 _jsi£ 
V sin 6»o 



d0- 



)(cos 2 + £) 



(B.10) 



Z COS 9n 



Vn 2 + z 2 



^nfsin- 1 



sm( 



sin 0o 



,sin 2 O , 



f2 sin 0q 

V^ 2 + z 2 



X = ^I(g,U/U ) 



(B.ll) 



where I(g,U/Uo) is defined in Eq. ( |3.1C| ). 

17 A check of our calculation is that the J = case reduces to Maldacena's calculation JlFj]. From Eq. (B.12), 
J = could mean two possibilities: 9q — > or l/fl — > oo. The former is correct: So — » gives the range 



AX = 7T- 



while I/O, — > oo gives Ax = 7rsecSo (valid only at So = 0). Therefore, the r2 — > limit 



is only a measure-zero subset of the So — > limit. We should take the So — * limit first to obtain J — > 0. 

Recalling that < < S o, the limit So — > ( J — > 0) makes the integrals on the left side of Eq. (BJ3) and 
the right side of Eq. ( |B.10|) equal, implying 



Ax 



X = l 



dy 



IF sin 



2 Vir/p,, V(f 2 -i)(y 2 + i-' 2 ) 

We defined the constant / in Eq. (B.4) such that it becomes Maldacena's £ Jl?! in the limit J — > 0. 
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The angular momentum and the regularized energy are given by 



j = — — / de 

k Jo 

V\ n 2 



sin t) cos I 



(sin 2 9 - sin 2 8) (cos 2 9 + 



(B.12) 



7T V^ 2 + / 2 



A" 



$7 sin # ( 



o 



Vn 2 + 1 2 



E 



$7 sin 6q 

VO 2 + 1 2 



7T 



(B.13) 



Parameters are determined by the equations matching the boundary conditions : 



0o 



dQ 



cos ( 







1 



sirr 9q 

which can be written as 

$7 sin 9q 



)(cos 2 8 + 



yr^T 2 " 



K 



V57 2 + 1 2 Wn 2 + i 2 



$7 sin 9, 



o 



dy 



y / iy 2_ 1){ y2 +g 2 ) 



I 2 K(ig) , 



(B.14) 



Ax 
2 



I 



n sin 6»o Jo 

I COS #0 

Vn 2 + 1 2 



d9- 



COS0a/(1 



sin z 



n ( sin 2 i 



sin 2 0q 

$7 sin 6*o \ 

Vn 2 + i 2 ) 



)(cos J 



+ IF 



L 
2~ 



I{9) ■ 



(B.15) 



(B.16) 



For given J and A%, #o an d ^ are determined from Eq (B.12) and ( |B.15| ). For these 
values, $7, I, and g are determined from Eq ( |B.14| ) and ( p3.6| ) . For this value of g, Uq is 
determined from Eq (B.16). Finally, u> is determined from the relation $7 = 

The right hand side of Eq. ( |B.14 ) is bounded from above and the left-hand side can be 
written as 18 ^J/sin0 o + it JcT^ cos Oy/l - sin 2 9/ sin 2 9 /^cos 2 8 + l 2 /Q 2 , where the 
integral in the second term is finite. Therefore the value of J is finite for any value of g and 
I. Since the string configuration in the X—U plane and the corresponding energy have the 
same expressions as in the Ax = case, J is finite when Uq touches the horizon and the 
interaction energy vanishes as g —* 0. In particular, at g = the maximum value 9* of 8q 

satisfies 

Vl — I 2 sinft,, \ 7r 

2 ' 



sin ( 



K 



y/l + I 2 sin 2 i 



y/l + I 2 Sin 2 

where we used f2 2 sin 2 9* = 1 — I 2 for g = 0. This occurs for finite J*, for any I, with J* 
reaching zero when I = 1 and A% = tt. By numerical computation, we find the curve in the 
(J, Ax) plane where E = 0; this curve is shown in Figure ^. 

18 This is compatible with the fact that J — > as 8q — > 0. 
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J/Va 

0.05 0.1 0.15 0.2 

Figure 6: Curve in the (J, A\) plane where E = (g = 0). 

C. Elliptic Integrals 

We use the following definitions of the elliptic integrals. 

r<p 

F(<p,k) = / -= 
Jo VI 



de 



sm 



E(<p,k) = [ V dOVl - k 2 sin 2 9 
Jo 

H(<p, n, k) 



de 



o (1 - n sin 2 6)y/l - k 2 sin 2 6 
K(k)=F{^,k" 



TT 



E(k)=E(-,k 
U(n,k) = n(|,n,fe) 

In some references and math-computing programs, F((p,k 2 ), E((p,k 2 ), Tl(ip,n,k 2 ), K(k 2 ), 
E(k 2 ) and II(n, /c 2 ) represent the same integral definitions. 
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